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Rapidly oscillating potentials with a vanishing time average have been used for a long time to 
trap charged particles in source-free regions. It has been argued that the motion of a particle inside 
such a potential can be approximately described by a time independent effective potential, which 
does not depend upon the initial phase of the oscillating potential. However, here we show that 
the motion of a particle and its trapping condition significantly depend upon this initial phase for 
arbitrarily high frequencies of the potential's oscillation. We explain this novel phenomenon by 
showing that the motion of a particle is determined by the effective potential stated in the literature 
only if its initial conditions are transformed according to a transformation which we show to 
significantly depend on the potential's initial phase for arbitrarily high frequencies. We confirm our 
theoretical findings by numerical simulations. Further, we demonstrate that the found phenomenon 
offers new ways to manipulate the dynamics of particles which are trapped by rapidly oscillat- 
ing potentials. Finally, we propose a simple experiment to verify the theoretical findings of this work. 
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I. INTRODUCTION 

The theoretical study of the dynamics of a particle 
moving in a rapidly oscillating potential has been an ac- 
tive field of research for the past 50 years [l|, 0, 0, 0, 0, 
H 0, B 0, 03 ■ In 1951, Kapitza showed that rapidly 
oscillating potentials can have a stabilizing effect on dy- 
namical systems jlj, which he demonstrated utilizing an 
inverted pendulum whose suspension point was forced to 
oscillate vertically. For high enough frequencies and large 
enough amplitudes of this oscillation the upwards posi- 
tion of the pendulum becomes stable. Kapitza showed 
that the pendulum feels a time independent effective po- 
tential that has a local minimum at the pendulum's up- 
wards position and thus imparts stability upon it. This 
phenomenon has since been termed "dynamic stabiliza- 
tion". In 1958, Paul used the principle of dynamic stabi- 
lization to trap charged ions by alternating electric fields 
fill ]. I° ns trapped within a Paul trap have and continue 
to be used in the development of frequency standards, 
other fundamental measures [l2| and quantum informa- 
tion jl3| ]. 

In order to grasp all the possible applications that are 
offered by atoms which are trapped by rapidly oscillating 
potentials, a descriptive theory of their dynamics is in- 
dispensable. For the motion of ions in a Paul trap such a 
theory is given by Mathieu's theory, which provides the 
exact solution of the ion's classical equations of motion 
in terms of the tabulated Mathieu functions [ill, EH EH • 
However, this theory is restricted to potentials of har- 
monic shape which are oscillated sinusoidally only. Be- 
sides this, Mathieu's theory does not raise intuition re- 
garding the atom's dynamics hindering the discovery of 
new phenomena. For these purposes, Kapitza's theory 
of the effective potential has thus far been suitable, suf- 
ficiently explaining the stabilization of the inverted pen- 
dulum [l|. This theory has since been generalized to ar- 
bitrary rapidly oscillating potentials @,E3| and extended 



to quantum mechanics [1, 0, E3| • 

In this paper, we describe a different phenomenon, 
which illustrates that the motion of a particle trapped by 
a rapidly oscillating potential with a vanishing time aver- 
age significantly depends on the potential's initial phase 
for arbitrarily high frequencies of its oscillation. We dis- 
cover that this phenomenon cannot be explained by the 
effective potential theory which we find to be incomplete. 
Therefore, we derive the necessary completion in order to 
describe this phenomenon. 

Kapitza's effective potential theory is based on a sep- 
aration of the particle's motion into a slow and a fast 
changing part. For high frequencies of the potential's os- 
cillation this separation can be performed in such a way 
that the slow part is governed by a time independent po- 
tential while the fast part has a vanishing time average 
and a negligible amplitude. In Kapitza's theory it is con- 
cluded that the particle's real motion is approximately 
described by its slow part, such that the real motion can 
be obtained by replacing its time dependent equation of 
motion by the time independent equation of motion of 
the slow part. However, doing so requires replacing the 
particle's initial conditions by the corresponding initial 
conditions of the slow part. If these two sets of initial 
conditions significantly differ from each other, the re- 
placement of the initial conditions is a crucial point in 
the determination of the real motion. In the context of 
adiabatic elimination a significant difference between the 
initial conditions of the full problem and the model prob- 
lem is called an "initial slip" [3ESE1- 

We show in our problem that there is an initial slip 
even for arbitrarily high frequencies of the potential's os- 
cillation. As a result, the approximate determination of 
the particle's real motion must be performed in two inde- 
pendent steps: Firstly the real initial conditions have to 
be transformed and then the time independent equation 
of motion of the slow part must be solved for these trans- 
formed initial conditions. In Kapitza's theory the trans- 
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formation of the initial conditions is omitted. We explic- 
itly derive this transformation and show that it signifi- 
cantly depends on the oscillating potential's initial phase 
for arbitrarily high frequencies. This explains the found 
phenomenon that the particle's motion depends on the 
initial phase. 

The phenomenon is clearly manifested in the particle's 
stability region (i.e. the region in phase space that con- 
tains all initial conditions of the particle leading to subse- 
quent trapping by the oscillating potential) for which we 
derive an analytical formula. Furthermore, we show that 
the phenomenon offers the possibility to manipulate par- 
ticles, that are trapped by rapidly oscillating potentials, 
in a simple and controlled manner. For example, modifi- 
cation of the oscillating potential's phase at some point 
in time by a discrete value (inducing a phase hop) results 
in a change of the trapped particle's average energy. We 
explicitly calculate this change in energy and show that 
it is significant for arbitrarily high frequencies of the po- 
tential's oscillation. For the special case of a sinusoidally 
oscillated arbitrary potential we also show that - in the 
framework of classical physics - it is significant even for 
arbitrarily low energies of the trapped particle. 

We confirm all theoretical findings in this paper by nu- 
merical simulations. For this purpose, we consider two 
simple examples: A particle moving inside an oscillat- 
ing Gaussian-shaped potential and inside an oscillating 
parabola-shaped potential. We also propose a simple ex- 
periment to validate our findings. 

Throughout the paper we use one-dimensional formal- 
ism, but the results presented here can be directly applied 
to two and three dimensions in cases where the motion is 
separable 0] • For high dimensional oscillating potentials 
where the particle's motion in not separable and in par- 
ticular chaotic, new phenomena may occur that are not 
discussed here [l9| . 

II. PARTICLE MOTION IN A RAPIDLY 
OSCILLATING POTENTIAL 

In this section, we derive an approximate theory for 
the motion of a particle that is moving inside a rapidly 
oscillating potential. This is done in two steps: Firstly, 
we formulate a model problem with a time independent 
"effective" equation of motion and then we show that this 
model problem possesses initial slip and derive the nec- 
essary transformation of the real problem's initial condi- 
tions to the corresponding initial conditions of the model. 
With the obtained theory we then verify the stated phe- 
nomenon of the dependence of the particle's motion on 
the initial phase of the oscillating potential. Finally, we 
describe a clear manifestation of this phenomenon and 
its main consequences. 

Newton's equation for a particle with mass m that is 
moving inside a one-dimensional rapidly oscillating po- 
tential V (x, t) can be written as 

mi = -V'(x, t) = -Vq(x) - V[(x)f(wt+<p), (1) 



where Vo(x) describes the time average part of V(x,t) 
and Vi(x) its oscillating part and f{ujf) its time modula- 
tion function which has period 2n, unity amplitude and 
vanishing time average over one period. The frequency 
of the potential's oscillation is denoted by ui. ipis the ini- 
tial phase of the oscillating potential V(x, t). Derivatives 
with respect to coordinates are denoted by primes and 
with respect to time by dots. The average over one period 
of oscillation (period-average) of a 27r-periodic function 
g(r) is denoted by a bar: g = ^ J n g{r) dr. If the fre- 
quency bj is very large compared to the other time scales 
of the system, the time dependent equation of motion ([1]) 
can be replaced by a time independent effective equation 
of motion that yields the period-averaged motion of the 
particle This effective equation of motion can be 

obtained by separating the motion x(t) of the particle 
into a sum of a slow part X(t), which depends on the 
slow time t, and a fast part £(r), which depends on the 
fast time r = cut [l[: 

x(t)=X(t)+t(r). (2) 

X(t) and £(t) will be referred to as the slow and the fast 
motion, respectively. For large frequencies uj, the ampli- 
tude of the fast motion £ can be assumed to be small, 
since due to its inertia, the particle does not have the 
time to react to the force which is induced by the oscil- 
lating potential, before this force changes sign. Accord- 
ing to [10], £ can for large frequencies thus be expanded 
in powers of the inverse frequency: £(r) = Y^Li ^^( T )- 
Requiring that on the fast time scale r, the fast motion £ 
is periodic with a vanishing period-average, the particle's 
period-averaged motion x(t) is given by the slow motion, 
x(t) = X(t). Further requiring that the definition of £ 
leads to a time independent equation of motion for the 
slow motion X(t), uniquely defines £ and yields the ef- 
fective equation of motion for X(f). This calculation is 
presented in Ref. [lOj. Its results in the leading order in 
lo^ 1 are the explicit expression of £, 

£(T,X,<p) = / (2)T [/(r+^)], (3) 

and the effective equation of motion 

mX = -V: S (X) (4) 
with the effective potential 

V cS (X) = V (X) + -A- 2 V[(Xf^j T [f{r)]j. (5) 

The standardized integral / t [/(t)] is the anti-derivative 
of f{r) with a vanishing period-average [20]. For 2tt- 
periodic functions /(r) that have a vanishing period- 
average and whose Fourier expansion is given by /(r) = 
En # o fne mT , it is defined as f[f{r)] = £„ #0 &/»e <nT - 
It can be applied repeatedly. Its multiple application (J 
times) is denoted by J^ T [/(r)] and its evaluation at the 



3 



point t by [f[- ■•]] To . 

The effective equation of motion ^ defines a simpli- 
fied model of the original problem. We now show that it 
possesses initial slip. For that matter, we compare both 
the position x(t) with its period-average X(t) and the 
velocity x(t) with its period-average X(t) determined by 
the original and the model problem, respectively. (For 
a numerical example see Sec. Ill, Fig. [3]) We restrict 
ourselves to the physically interesting case of a particle 
that is trapped by a rapidly oscillating potential with a 
vanishing time average (i.e. Vq = 0). The position X(t) 
and the velocity X(t) of the period-averaged motion are 
determined by equation ([4J. Thus, for the particle to 
be trapped, the effective potential ((5|) must have a local 
minimum. The spatial width of this local minimum de- 
termines the spatial confinement of the particle, and is, 
according to Eq. J5]), independent of lo. Thus, the posi- 
tion X{t) of the trapped particle's slow motion is of the 
order co°. According to ((3]), the position £(i) of its fast 
motion is of the order to~ 2 and therefore in the limit of 
large to negligible compared to the position X(t) of its 
slow motion. This means that, for large lo, one can ap- 
proximate x(t) = X(t) + £(t) w X{t) (see Fig. da)). 

The depth of the effective potential's local minimum 
determines which kinetic energies the trapped particle 
can have. It is, according to Eq. ([5]), of the order u~ 2 . 
Thus, the velocity X(t) of the particle's slow motion is 
of the order lo" 1 . The velocity £(i) of its fast motion is 
given by the derivative of Eq. with respect to time, 



£(*,¥>) = - — V{{x{t)) 
mio 



lf(T + ¥>)} 



, (6) 



which is a term of the order w _1 and thus of the same 
order as X(t). Therefore, the velocity £(t) of the trapped 
particle's fast motion is for large lo not negligible com- 
pared to the velocity X(t) of its slow motion. As a result, 
even for large w, it is x(t) = X(t) +£(i) 76 X(t) (see Fig. 
[DJb)), such that in general it is i(0) 76 X(0). This shows 
that there is an initial slip. 

This observation explains that in order to determine 
the particle's period-averaged motion with the effective 
equation of motion ^ for given initial conditions x(0) 
and x(0) the initial condition transformation 



X(0) = x(0), 
X(Q) = i(0)-£(0), 



with 



m = v ; (xm 



(7) 



(8) 



T = 



has to be performed, whereas the contribution of £(0) in 
Eq. ([7]) cannot be neglected. 

Expression (|8|) depends on the term [J T [/(t+(/?)]] , 
which is a function of the initial phase (p of the oscil- 
lating potential V(x,t). The initial condition transfor- 
mation ([7]) therefore depends on ip, which implies that 



the resulting period-averaged motion X(t) of the parti- 
cle also depends on ip (For a numerical example, see Sec. 
Ill, Fig. 13]). Further inspection of expression ([8]) shows, 
that it also depends on the initial position x(Q) of the 
particle. This implies that the degree of dependence of 
the particle's motion on the initial phase is determined by 
the particle's initial position. In the special case that the 
term [/ T [/( T + ( i 5 )]] T _ vanishes, the initial conditions of 
the slow motion, X(0),X(0), equal the real initial condi- 
tions, x(0), i(0), and there is no initial slip. In this case, 
the motion of the particle is determined by the time in- 
dependent effective potential © alone. Since / t [/(t)] is 
27r-periodic and has vanishing average, there always exist 
two phases pi,(p2 £ [0, 2ir] for which [J [f(r+ipi)]\ T 
and [J t [/(t+952)]] r _ vanish. Thus, for every given 
rapidly oscillating potential, there exist two phases for 
which the motion of any particle inside this potential is 
determined only by its time-independent effective poten- 
tial. 

The significance of the highlighted dependence of a 
particle's motion on the initial phase is manifested in the 
particle's stability region, which illustrates the particle's 
trapping condition. The stability region is the smallest 
region in phase space that contains all initial conditions 
of the particle leading to subsequent trapping by the os- 
cillating potential. Since the effective potential ([5]) in 
equation (U) is time independent, the principle of energy 
conservation can be invoked for its calculation. For a 
trapping effective potential we define X = as the po- 
sition of the bottom of the trap and X = L as the po- 
sition of the trap boundary. Thus, the initial conditions 
X(0), X(0) belong to the stability region if they obey the 
inequality 



-roX(O) 2 + V eS (X(0)) < V cS (L). 



(9) 



Substituting Eq. (UJ, solving for x(0), and neglecting 
terms of higher order than uj^ 1 yields 



i(0) ^ 



-y((x(o)) 



\f(r+<p)] 



± 



^([[f{T)]j^jv{{Ly-v{{xm\ (10) 

which defines the stability region for large lo. The right- 
hand side of Eq. (TTUJ) consists of a phase dependent and 
a phase independent term which are both of the same or- 
der, such that the stability region significantly depends 
on the phase if for arbitrarily large lo. In Sec. Ill, we 
explicitly calculate the stability region for different ip for 
two numerical examples. 

In the preceding paragraphs we showed that, due to 
the dependence of the fast motion's initial velocity £(0) 
on the initial phase tp and due to the significant coupling 
between the particle's slow motion X(t) and £(0), the 
particle's slow motion X(t) is significantly coupled to <p. 
ip can therefore be used to manipulate the particle's mo- 
tion. For example, if tp is modified instantaneously by a 
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discrete value Aip (a phase hop), the position X(t) of the 
particle's slow motion will, on the time scale of the slow 
motion, change in a non-smooth manner and the velocity 
X(t) will, on this time scale, change non-continuously by 
a discrete value, AX. According to Eq. (JT]), this change 
is given by 



Ai(t P h) 



1 



Vi(X(t ph )) ■ A 



with 



A 



[f(T+<p+A<p)] 



[/(r+v)] 



(11) 



. (12) 



where i p h denotes the time at which the phase hop is 
induced. This change in the velocity of the particle's 
slow motion results in an instantaneous change AE (on 
the slow time scale) of the period-averaged total energy 
E of the particle, which is given by 

AE(t ph ) = mX(t ph )AX(t ph ) + ^m(AX{t ph )) 2 , (13) 

where A(£ p h) denotes the velocity of the slow motion 
right before the phase hop. As can be seen from equa- 
tions (fTTj) and (fT3"|) , AE(t p b) is of the same order (w~ 2 ) 
as E. Thus, the period-averaged total energy E of a 
trapped particle can be significantly manipulated by a 
phase hop for arbitrarily large u>. 

As an example, we calculate AE for the special case 
f(r + (p) = cos(r + 9?), and assume that the phase hop 
is induced at a time t = i p h for which J T [cos(r + 
<y5)]r=wt ph = (which is the case twice each period 
of the fast motion), and for which the kinetic energy 
of the slow motion equals three-quarters of the total 
period-averaged energy, that is ^mXftph) 2 = jE (this 
is the case four times each period of the slow motion). 
Due to energy conservation of the slow motion, it is 
V cS (X(t ph )) = ^V({X{t ph )) 2 = \E. We further as- 
sume that the size of the phase hop A<p equals ? , such 
that f T [cos(t + (p + A(p)] T=LO t ph = 1. Then 



AE(t ph ) 



mX(t ph ) — Vl(X(t ph )) + — ^ 




V[{X{t ph )f 
(14) 



where the minus sign holds for X(t p h)V{(X(t p h)) < 
and the plus sign for X(t ph )V{(X(t ph )) > 0. AE is 
therefore independent of u) and proportional to E, so that 
the relative change, =^-, is independent of E (note, that 
this holds only in the limit of large oj and small E) . It is 
therefore possible to take away more than two thirds of a 
trapped particle's energy or to give more than five thirds 
of its energy by inducing a phase hop for arbitrarily high 
frequencies and - in the framework of classical physics - 
for arbitrarily low energies (for numerical examples see 
Sec. Ill, Fig. 2] and Fig. 0. In Sec. IV, we propose an 
experiment which may verify this result. 



III. NUMERICAL EXAMPLES 

In this section, we consider two simple examples in 
order to demonstrate the phenomenon found in the pre- 
ceding section and to compare the derived analytical pre- 
dictions to results of numerical simulations. The first 
example to be considered is a particle moving inside an 
oscillating Gaussian-shaped potential with a vanishing 
time average of the form 



2x 2 " 

V u (x, t) — 7exp ( J" ) cos(wi + ip), 

w 



(15) 



which is illustrated in Fig. [Ha). Potential (TT5"|) defines 
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FIG. 1: (Color online) (a) Oscillating Gaussian-shaped poten- 
tial (|15|l and (b) oscillating clipped parabola-shaped potential 



two time scales: The time scale associated with the fre- 
quency ui of the potential's oscillation and the time scale 
associated with the oscillation frequency lo osc of a low 
energy particle in the minimum instantaneous potential. 
The theory which was derived in Sec. II is valid for large 
frequencies to compared to the frequencies which are as- 
sociated with the other time scales of the system. In the 
following, we therefore consider the regime w>w osc . 

In the preceding section, we showed theoretically, that 
for correctly chosen initial conditions the position of a 
trapped particle can be approximated by the position of 
its slow motion, whereas its velocity cannot be approxi- 
mated by the velocity of its slow motion. To verify this, 
we numerically simulate a particle's trajectory in position 
and velocity space and compare it to its period-average. 
This is done by numerically integrating the time depen- 
dent Newton equation 



47 

mx = — ,a;exp 
"0 



2x 2 " 

— 2" ) cos(wi + if). 



(16) 



The results are shown in Fig. [5] In the figure, one can 
clearly see that the trajectory in position space is well 
described by its period- average (Fig. Efa)), whereas in 
velocity space (Fig.[2jb)) this is not the case. Additional 
numerical simulations showed that this result does not 
change when ui is increased, in agreement with the de- 
rived theory. 

In Sec. II, we showed that the significant difference be- 
tween the particle's velocity and its period-average (Fig. 
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result is shown in Fig. [3] The charted trajectories appre- 
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FIG. 2: (Color online) Trajectory and period-averaged tra- 
jectory of a particle inside the oscillating Gaussian-shaped 
potential (|15(l obtained from numerical simulations (a) in po- 
sition space and (b) in velocity space. Solid curve (red): 
real trajectory, dashed curve (black): period-averaged tra- 
jectory. Parameters: u) — 7tJ sc; <p — 0. Initial conditions: 
x(0) = 0.32 Wo', x(0) — 0. The real and period-averaged tra- 
jectory are almost identical in position space but largely differ 
in velocity space. 



[2Kb) ) leads to a dependence of the particle's motion on 
the initial phase tp of the oscillating potential. To ver- 
ify this conclusion, we numerically simulate the particle's 
trajectory in position space for different ip for a given fre- 
quency u) and given initial conditions of the particle. The 
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FIG. 3: (Color online) Trajectories in position space of a par- 
ticle inside the oscillating Gaussian-shaped potential (|15|) for 
different ip, obtained from numerical simulations. Parameters: 
lu = 7u sc and ip — (red curve), ip = 2 (blue curve), tp = n 
(also red curve), and <£ = %- (purple curve). Initial conditions: 
x(0) = 0.32 wo; x(0) — 0. The shown trajectories are almost 
indistinguishable from their theoretical predictions calculated 
from Eqs. ([5]) and Q, which are therefore not shown in the 
figure. 

ciably differ from each other, as predicted by the derived 
theory. On the slow time scale they have different phases 
and different amplitudes. Figure [3] shows as well, that all 
trajectories appear to be generated from the same time 
independent effective potential. Since the trajectories for 
different tp are not the same although their initial condi- 
tions are, this verifies that the trajectories are determined 
in addition by a transformation of their initial conditions 
which depends on ip. The shown trajectories are almost 
indistinguishable from their theoretical predictions cal- 
culated from equations (O and ((7|), which are therefore 
not shown in the figure. 

Another result of Sec. II to be verified by numeri- 
cal simulations is the prediction of the particle's motion 
when phase hops in the oscillating potential's time mod- 
ulation function are induced. Figure 2] shows the nu- 
merically simulated trajectory of the particle for a phase 
hop of size Ap> = J induced at a time t — £ p h for 
which J T [cos(r + ^)] r=aJtph = 0, ±mX(i ph ) 2 « f and 

X(t p h) V{(X{tph)) > 0. In agreement with the analytical 
prediction (fH)) . the period-averaged energy of the parti- 
cle in Fig. [4] indeed changes by a factor of approximately 
2.7. Furthermore, the shown trajectory is almost indis- 
tinguishable from its theoretical prediction (not shown in 
the figure) calculated from equations (0 and (J7J), whereas 
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FIG. 4: (Color online) Trajectory in position space of a parti- 
cle inside the oscillating Gaussian-shaped potential (|15[1 when 
a phase hop is induced, obtained from numerical simulations. 
Solid curve (red): real trajectory, dotted curve (black): imag- 
inary continuation of trajectory when no phase hop was in- 
duced. Parameters: u == 7tJ sc, <p = 7T, time of phase hop 
tph = 96 • — , and phase hop size A(p ===== -|. Initial condi- 
tions: x(0) = 0.32 wo; i(0) = 0. The phase hop changes 
the particle's total energy by a factor of approximately 2.7 
in agreement with equation (I14|l . The shown trajectory is 
almost indistinguishable from its theoretical prediction (not 
shown) calculated from Eqs. (JS| and Q. 



the part of the trajectory after the phase hop is obtained 
from these by assigning the time of the phase hop to the 
initial time. 

Figure O shows the numerically simulated trajectory 
of the particle when several phase hops of the size 

U = 



are induced. The moments of the phase hops 



were chosen such that J T [cos(t + (p)] T _ 



0, \mX{t%f « and X(t%>)V{(X{t%)) 



"ph 
< 



0, 



whereas X(t^l), X(t^l), and E(t^l) denote the position, 
the velocity, and the energy of the particle's slow motion 
right before the ith phase hop, respectively. In agreement 
with the analytical prediction (fr4|) . the period-averaged 
energy of the particle in Fig. [5] changes in total by a fac- 
tor of approximately j^o ■ Figure [5] demonstrates that 
this substantial change in energy can be achieved even 
when the phase hops are induced very quickly one after 
another (within a fraction of the period of the particle's 
slow motion). The shown trajectory is almost indistin- 
guishable from its theoretical prediction (not shown) cal- 
culated from Eqs. (J5]) and ([7]). 

The last theoretical prediction of Sec. II that is to be 
verified by numerical simulations in the following is the 
derived equation (|10p for the stability region. For the 
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FIG. 5: (Color online) Trajectory in position space of a parti- 
cle inside the oscillating Gaussian-shaped potential (|15[l when 
a series of four phase hops is induced, obtained from numerical 
simulations. Solid curve (red): real trajectory, dotted curves 
(black) : imaginary continuations of trajectory when the phase 
hops were not induced. Parameters: cu = 7 ui osc , <P = k, times 



of phase hops t 



(4) 

ph 



(i) . 
ph 



= 84 



27T ,(2) . 
w ' ph ' 



= 92.75 



$ = 103.5- 



112.25 • ^j, and phase hop size Aip == ?. Initial con- 
ditions: x(0) == 0.32 wo; x(0) = 0. The phase hops change 
the particle's total energy by a factor of approximately 1/170 
in agreement with equation (|14fl . The shown trajectory is 
almost indistinguishable from its theoretical prediction (not 
shown) calculated from Eqs. ((5| and lff[l. 



particle inside the oscillating Gaussian-shaped potential 
(fT5f . it yields 



1 4 7 



-jxexp 



simp 



1 V87 



^exp(-l)-x2exp (17) 



A graph of this analytical prediction of the stability re- 
gion for different tp is shown in Fig. [S] (solid curves) . Be- 
fore comparing it to the results of numerical simulations, 
we first discuss it separately. In Fig.[6l one can clearly see 
that the analytically predicted stability regions for differ- 
ent ip significantly differ from each other. For the initial 
phases ip = and tp = Tr, the stability regions are identical 
and axially symmetrical to the coordinate axes. For the 
phases y = f and cp = =y , the stability regions are not ax- 
ially symmetrical to the coordinate axes and particularly 
not to the position axis. This demonstrates, that the 
motion of the particle cannot be determined by a time 
independent effective potential alone, since the stability 
region of any time independent potential is due to en- 
ergy conservation axially symmetric to the position axis. 
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FIG. 6: (Color online) Stability regions for a particle inside 
the oscillating Gaussian-shaped potential (|15p for different 
phases ip and for uj = 1 uj sc- Solid curves: Analytical pre- 
dictions ([17)1 for (p — (red), <p=\ (blue), (p = n (also red), 
and ip = (purple). Dashed curves: Results of numerical 



simulations for tp = (red), <p = 
Dotted curve: For <p = ir (red). 



(blue) , and ip = ^ (purple) 



Compared to the stability regions of the phases p = 
and ip = Tr, those of the phases ip = | and tp = ^ further 
look distorted and rotated around the origin. One can 
see, that the area of the stability region is independent 
of the initial phase p, which immediately emerges from 
equation (fTO]) . The overlap of two stability regions of dif- 
ferent ip in units of their total area, however, depends on 
the initial phase of the respective overlapping stability 
regions, but is according to Eq. (I10[) independent of w. 
The spatial width of the stability region is independent 
of the initial phase ip and according to Eq. (fTU)) also inde- 
pendent of cu. Finally, Fig. [B] shows that the boundaries 
of all stability regions intersect at x = 0, and thus the 
escape velocity of the particle that is initially situated at 
x = is independent of the initial phase w, in agreement 
with the numerical results presented in [21] . This is, be- 
cause for x(0) = the initial condition transformation ([7]) 
is independent of <p, since it is V((0) — Q. 

We now compare the analytical prediction for the sta- 
bility region (|17p to results of numerical simulations, 
which are shown as dashed and dotted curves in Fig. 
IH1 The numerical simulations were performed as follows: 
First, the particle's phase space was discretized into a 
grid of equally spaced points, each of which define pos- 
sible initial conditions of the particle. Then, the time 
dependent Newton equation (|16p was integrated numeri- 
cally for each of these points. Thus, all the points, which 
generated a trajectory that stayed within a prespecified 



spatial region within a prespecified large time, were de- 
fined as points of the stability region. The prespecified 
values for the spatial region (|x| < ^) and the large time 
(t w 200 • ^j) were chosen such that further increase in 
these values did not change the result appreciably. The 
step size was chosen to be At = t^k ■ — , since further 

" 100 Ul 1 

decrease of the step size did not alter the results appre- 
ciably. Figure [6] shows a very good agreement between 
theory and simulations. Additional simulations showed, 
that this agreement increases when w is increased, as ex- 
pected from the theory. 

So far we have considered a particle moving inside the 
oscillating Gaussian-shaped potential (fT5|) . Another ex- 
ample to be considered in the following is a particle mov- 
ing inside an oscillating clipped parabola-shaped poten- 
tial with a vanishing time average of the form 



V cp (x,t) = Fx (2) cos(wt + </>), 



with 



Vi(x) 



-kx 2 + kL 2 for \x\ < L, 
else, 



(18) 



(19) 



which is illustrated in Fig. [ljb) . This is a problem able to 
be solved exactly, because the particle's time dependent 
Newton equation 



mx 



2kx cos(o>i + ip) 




for \x\ < L, 
else 



(20) 



is for I a; I < L equivalent to the homogeneous Mathieu 
equation, whose solution i s g iven by a linear combination 
of the Mathieu functions [lj, EH Gil > an( i f° r it is 

equivalent to the trivial differential equation x = of free 
motion. Therefore, the considered example allows for a 
comparison between the theory that was derived in Sec. II 
and the Mathieu theory. We choose L = ^ and k 

so that in Eq. (I18[) the reference frequency uj osc ■ 



. 7 
— 



as well as the maximum instantaneous potential depth 7 
are the same as in (fl~5|) . In the following, we calculate 
the particle's stability region - once by using the derived 
equation p0[) and once by applying Mathieu's theory. 
Equation (flT)]) yields 



!.(cp) 



5 



1 - 

— |i sin(y) ± 



1 2\/2< 



(21) 



A graph of this analytical prediction of the stability re- 
gion is shown in Fig. [7] (solid curves) for different ip. 
Figure [7] also shows the results obtained from Mathieu's 
theory (dashed and dotted curves). It can be seen that 
the stability regions of the particle have the shape of an 
ellipse. Figure [7] shows a very good agreement of our 
derived theory and the Mathieu theory. In order to con- 
firm the accuracy of the numerical simulations which we 
performed for the oscillating Gaussian-shaped potential 
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FIG. 7: (Color online) Stability regions for a particle inside 
the oscillating clipped parabola-shaped potential (|18p for dif- 
ferent phases <p and for uj — 7 u sc- Solid curves: Analytical 
predictions (|2ip for p = (red), <p — 5 (blue), tp = n (also red), 
and ip = (purple). Dashed curves: Results of Mathieu's 
theory for <p = (red), V 9 = § (blue), and <p — 4^ (purple). 
Dotted curve: For <,9 = 7r (red). 



in the preceding paragraph, we repeated these simula- 
tions for the parabola-shaped potential finding that the 
obtained results were identical to those obtained from 
Mathieu's theory. All main results which we demon- 
strated for the oscillating Gaussian-shaped potential were 
found also for the oscillating parabola-shaped potential. 



IV. PROPOSAL OF AN EXPERIMENT 

In this section, we propose an experiment to verify our 
theoretical results. Consider an ensemble of atoms that is 
trapped inside a rapidly oscillating potential with a van- 
ishing time average, whose amplitude is chosen such that 
the ensemble uniformly fills its stability region. Hence 
the fraction of atoms which remain trapped after a phase 
hop of the oscillating potential is induced is equal to the 
mutual overlap of the stability region with the initial time 
given by the time right before the phase hop and the sta- 
bility region with the initial time given by the time right 
after the phase hop. This overlap is exactly determined 
by equation (|10[) . It is a function of the time i p h when 
the phase hop is induced and of the size Atp of the phase 
hop. Figure [5] shows this overlap for a two-dimensional 
oscillating Gaussian-shaped potential. One can see that, 
regardless of i p h, at least 70% of the atoms will leave the 
trap for appropriate Atp. Therefore, when measuring the 
ratio between the numbers of trapped atoms before and 




0.00 -h 1 r- 1 

7t 2 71 

Acp [rad] 

FIG. 8: (Color online) Overlap of the stability region of ini- 
tial phase ip with the stability region of initial phase ip + Aip 
for different t p h as a function of Aip obtained from the 
two-dimensional extension of equation (|10[l for an oscillat- 
ing Gaussian-shaped potential. Solid curve (red): (u}t p \ 1 + 
ip) mod 2tt—^, dashed curve (purple): (ait^-Hp) mod 27r = 7r, 
dotted curve (blue): (cotp^ + p) mod 2tv—^-. The overlap of 
the stability regions is given in units of their total volume. 



after the phase hop for different phase hop sizes, one can 
expect a clear signal from the experiment. 



V. CONCLUSION 

In conclusion, we investigated the classical dynamics 
of a particle moving inside a rapidly oscillating potential 
with a vanishing average. We found, that the motion 
of such a particle significantly depends on the oscillating 
potential's initial phase in the case that the particle is 
trapped. The effective potential theory described in the 
literature failed to describe this phenomenon. We found 
that this theory is incomplete. It states that the motion 
of the particle is determined by a time independent ef- 
fective potential, whereas we could show that this is true 
only if the particle's initial conditions are transformed ac- 
cording to a non-negligible transformation which is sig- 
nificant even for arbitrarily high frequencies of the po- 
tential's oscillation. We derived this transformation and 
showed that it depends on the potential's initial phase 
which explains the found phenomenon. We also showed 
that there always exist two initial phases of the oscillating 
potential such that the transformation can be neglected. 
Therefore, the effective potential theory in its previous 
form, which omits this transformation, can only describe 
these two special cases correctly. The phase dependence 
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offers a new possibility to significantly manipulate the 
dynamics of a trapped particle. We explicitly calculated 
the change in a particle's energy caused by a phase hop 
in the potential's modulation function. We found that 
it is significant for arbitrarily high frequencies of the po- 
tential's oscillation and - in the framework of classical 
physics - for arbitrarily low energies of the particle. It is 
the subject of future research to find applications of this 
novel tool of particle manipulation. The results presented 
in this article, should also be extended to the quantum 
regime. 
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